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FOREWORD 


Soon after the Mexican volcano El Chichdn erupted during the period March 28 to April 4, 
1982— spewing large amounts of ash and gaseous products into the atmosphere— unusual effects 
began to appear in certain geophysical parameters determined from satellites using remote sensing 
techniques. For example, sea surface temperatures obtained in equatorial regions from the NOAA-7 
operational satellite were observed to be too low by several degrees (°C) when compared to ship and 
buoy measurements, and total ozone retrievals from the Nimbus-7 satellite were observed to be 
unreasonably large in the vicinity of the volcanic cloud. 

As a result, an ad hoc El Chichdn Remote Sensing Working Group (ERSWG) was formed, 
consisting of scientists and engineers associated with some currently active satellite sensors, to 
exchange information on how the eruption was affecting the various sensors and to coordinate 
activities with respect to developing a more complete understanding of the event and to formulating 
additional plans for further research. 

The purpose of this document is to combine under one cover a number of individual reports 
concerning the history of F.1 Chichdn eruption, a treatment of some relevant chemical and radi- 
ative processes involved, and a discussion of the expected and observed effects of the volcanic cloud 
on 10 currently active satellite sensors. 

Coming relatively soon after the eruption, many of the discussions in this document are indeed 
(as the title implies) preliminary and incomplete. However, it was deemed advisable to disseminate 
this information as soon as possible in order to inform other workers in the field of these activities 
and to provide a point of departure for further research that will undoubtedly continue for a long 
time to come on the effects of volcanic eruptions on remote sensing. 

William R. Bandeen 
Chairman, ERSWG 
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ABSTRACT 


The aerosols and gases resulting from the El Chichdn volcanic eruption had. and may still have, 
significant effects on satellite measurements of the properties of the earth’s atmosphere. The sea 
surface temperature measured by the AVHRR was biased up to -2.5°C for many months. The total 
amount of ozone derived from TOMS with a standard algorithm was greatly in excess of the possible 
value. This apparent excess can now be explained in terms of additional absorption by S0 2 - 
Infrared temperature sounders have observed both positive and negative anomalies. These effects 
and others on many satellite measurements are addressed following a discussion of the history and 
composition of the ejecta remaining in the stratosphere. Finally, recommendations are made for 
further study to account for the effects of volcanic eruptions on satellite observations and for use 
of such observations to measure the characteristics of the ejecta. 
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RADIATIVE EFFECTS OF THE 
EL CHICHON VOLCANIC ERUPTION: 

PRELIMINARY RESULTS CONCERNING REMOTE SENSING 

1. INTRODUCTION 

Satellite measurements of properties of the earth and its atmosphere were seriously affected 
by the debris thrown into the atmosphere as a result of the massive El Chichon volcanic eruptions 
during the period March 28 to April 4, 1982. The volcano is located in Mexico at latitude 17.3°N 
and longitude 93. 2° W (Fig. 1.1). Spurious estimates of total atmospheric ozone from satellite 
measurements were soon apparent, and difficulty in retrieving sea surface temperatures from in- 
frared measurements has persisted. Many other satellite observations are affected, also. The purpose 
of this report is to make a preliminary assessment of the effects of the El Chichdn atmospheric 
contamination on most satellite current observations of the earth and its atmosphere. 

The volcanic aerosols residing in the stratosphere change the climate by modifying the radiation 
balance of the earth-atmosphere system. Relevant parameters resulting from the recent Mount 
St. Helens’ eruptions have been extensively documented by Deepak (1982), in Science, v. 21 1, 
February, 1981, and by Newell and Deepak (1982). and for the Soufriere eruptions in Science, 
v. 216, June 4, 1982. The El Chichon aerosol is expected to cause a stronger climate effect, because 
more aerosols occur higher in the stratosphere where they will last longer. The current satellite 
observations will measure some characteristics of this stratospheric layer on a global scale. 

This rerort is organized to follow with a discussion of the history of recent El Chichon erup- 
tions, largely based on satellite observations. The important chemical processes in the stratospheric 
volcanic cloud will be reviewed. Their radiative properties and also those of the aerosols persisting 
in the stratosphere are discussed. Then a preliminary assessment of the effect of the El Chichon 
debris is made for the following satellite sensors: 


Mi 
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Advanced Very High Resolution Radiometer (AVHRR) on NOAA-6 and 7, TIROS Operational Vertical 
Sounder (TOVS) on TIROS-N, Visible and Infrared Spin Scan Radiometer (V1SSR) and the V1SSR 
Atmospheric Sounder ( VAS) on the Geostationary Operational Environmental Satellite (GOES), 
Solar Backscattered Ultraviolet/Total Ozone Mapping System (SBUV/TOMS) on Nimbus-7, Earth 
Radiation Budget (ERB) on Nimbus-7, .oastal Zone Color Scanner (CZCS) on Nimbus-7, Scanning 
Multichannel Microwave Radiometer (SMMR; on Nimbus-7, and the Multispectral Scanner (MSS) 
and Thematic Mapper (TM) on Landsat. The authors of each subset tion appear it its end. Con- 
clusions and recommendations to reduce errors in the derived parameters or for deriving new 
information concerning the volcanic debris are given. 

References (Section 1) 

Deepak, A. (Editor), 1982: Atmospheric Effects and Potential Climate Impact of the 1980 Erup- 
tions of Mount St. Helens. NASA Conference Publication 2240, Scientific and Technical 
Information Branch, 303 pp. 

Newell. R. E. and A. Deepak, 1982: Mount St. Helens Eruptions of 1980. NASA, Scientific and 
Technical Information Branch, Washington, D.C. 1 19 pp. 
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2. ERUPTIONS OF EL CHICHON VOLCANO 

2 i. GOES AND NOAA-6 OBSERVATIONS OF EL CHICHON ERUPTIONS 

On March 28, 1982, at 2332 local time, (add 6 hours to change to GMT) El Chichon 
volcano erupted violently, the first of four major eruptions. The GOES first recorded the eruption 
at midnight Imagery from this satellite is acquired approximately every 30 minutes in both the 
T isible (0.55-0.75 micrometers) and thermal infrared (10.5*12.5 micrometers) at a maximum 
. solution of 1 and 8 km, respectively. On March 29, at 0230, satellite digital data gave 
lowest eruption cloud brightness temperatures of -75.2C. Comparison with radiosonde data 
at Veracruz, Mexico (19.15N, 96.12W), showed that this temperature was readied at a height of 
1 6.3 km. The Veracruz data also showed that the altitude of the tropopause was 16.5 km, thus 
the eruption column probably did reach the tropopause. 

By 0800 the ash cloud began to dissipate, the minimum satellite-derived brightness temper- 
atures being -63.2C at this time. The temperature corresponded to a height of 13.3 km. The ash 
cloud spread m a northeasterly direction, although a westerly component over the Pacific was 
a lso evident (Fig. 2.1). Upper air wind data at Veracruz showed winds from the SW between 10.4 
and 1 5.2 km, and winds from the NE at 24 km. Thus the easterly drift was tropospheric and the 
Pacific com.onent was stratospheric. By the end of the day the main ash cloud stretched from 
Cuba to 103Wand covred an area of approximately 600,000 km"\ 

A f °2 0 on April 3, 1982, El Chichon erupted again and the first GOES satellite recording of 
the „</ent occurred at 0300. This eruption was not as large as the March 28 event Digital data at 
0400 showed minimum brightness temperatures of -7 1 .2C, corresponding to an altitude of 1 5. 1 km 
when compared to Veracruz radiosonde data. This was below the tropopause which was at 17.0 km. 
At 0800 - -5 9. 2C was the lowest brightness temperature recorded by the GOES data, corresponding to 
a h-ight of 1 2.5 km. After the eruption the ash cloud spread both NE and SW (Fig. 2.2) at heights 
of between 13.7 to 16.3 km and 0 to 4.9 km, respectively, compared with Veracruz upper air data. 

2-1 
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ASH CLOUDI 


^STRATOSPHERIC CLOUD 


Figure 2. 1 . C'iOLS enhanced thermal IR image of El Chichon ash cloud on 
March 2d. 1^82 at I 500 local time. The cloud is from the 
March 28 eruption at 2332 local time. 
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Figure 2.2 GOFS unenhaneed thermal IR image of FI Chichon ash cloud, which was in the troposphere, 
on April 3. 1982 at 0900 local time. The cloud is from the second eruption at 0-50 local time. 




Late on April 3, 1982, and again early on April 4, 1982, there were two additional major 
eruptions of El Chichon, the second occurring within 10 hours of the first. The first eruption began 
about 2000 and the subsequent ash cloud spread over southern Mexico, Belize, and northern 
juatemala (Fig. 2.3), elongated NE-SW. Upper air data at Veracruz showed stratospheric northeast 
winds from 20.7 to 24.0 km and tropospheric southwest winds from 13.7 to 18.5 km. At 2200 
minimum satellite brightness temperatures of -78. 2C were recorded, lower than the lowest radiosonde 
temperature of -77.5C at the tropopause (16.9 km). On April 4 at 0030 minimum satellite bright- 
ness temperatures were -71.2C, corresponding to an approximate height of 15 km. 

The final and largest eruption of El Chichon began at 0522 on April 4 (Fig. 2.4). The NOAA-6 
polar-orbiting satellite recorded the eruption at 0744 in both the visible band (0.58-0.68 pm) 
and thermal IR band (10.5-1 1.5 micrometers) at a -esolution of 1 km. The thermal IR data shewed 
satellite-brightness temperatures at least as low as -83.0C, a temperature below the lowest Veracruz 
radiosonde temperature of-75.3C at the tropopause (1 6.2 km). This eruption clearly penetrated the 
tropopause and was primarily responsible for the stratospheric cloud that subsequently circled t‘.e 
Earth; the Beginnings of which are visible drifting westward in Fig. 2.5. 

M. Matson 
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Figure 2.3. GOES unenhanced thermal IR image of El Chichon ash cloud between two eruptions on 
April 4, 1982 at 0100 local time. The cloud is from the third eruption on April 3 at 2000 local time. 






Figure 2.4. GOES unenhanceu thermal IR image of the fourth El Chichon eruption on April 4. 1982 
at 0600 local time. The ash cloud is from the third eruption on April 3. 
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Figure 2 5. GOES unenhanced thermal I R image of El Chichon ash cloud on April 4, 1982 

at 1730 local time, after the final eruption. 


2.2. OBSERVATIONS OF EL CHICHON CLOUD CHARACTERISTICS 

The eruption of El Chichon is moderately well documented compared to earlier major erup- 
tions such as Agung and Fuego because of its relative proximity to the U.S., the availability of 
advanced satellite sensors, and a fortuitous passage of the early stratospheric cloud over the NOAA 
Mauna Loa station. In this section the various observing platforms and instruments are indicated 
and the initial findings are described. The SEAN Bulletins (Smithsonian Institution, 1982) contain 
additional information about the eruptions. 

2.2.1. Observing capabilities 

2.2. 1.1. Operational satellites 

The ash and liquid phase clouds were initially sighted from the GOES and NOAA-6 
and 7 spacecraft using visible and infrared channels. The visible light observations provided 
information about the drift and dispersion of the clouds as described in the previous section. 
In addition, estimates of plume top height were derived from IR window channel temper- 
atures as described in Section 2.1. 

Effects of the cloud were found in sea surface temperatures derived from AVHRR data. 
Anomalously low values relative to s’uip observations were used to map the drift and spread of the 
cloud. Fig. 4.2a illustrates the affecttd areas during the month of May. These results are con- 
sidered in more detail in section 4.1. 

2. 2. 1.2. Research satellites 

The primary research satellites in operation during El Chichon eruption were Nimbus 7 
and Solar Mesospheric Explorer (SME). The sensors on Nimbus 7 of interest for eruption studies 
are SBUV/TOMS, Stratospheric and Mesospheric Sounder (SAMS), SMMR, ERB. and CZCS. The 
effects of the cloud as a contaminant to the radiances observed with the Nimbus sensors, except 
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for SAMS, are considered in detail in section 4. Nimbus data processing is generally delayed until 
definitive instrument parameters and validation studies are complete. However, in the case of 
TOMS the data were specially processed in near real time during the eruptions. Thus it was 
possible to observe the signature of the cloud immediately and to track its movement for guidance 
in sampling aircraft flight planning. 

Strong absorption was found in the two shortest wavelength channels of 1 (Krueger, 
1983). This absorption, which mimicked that of ozone in the retrievals, was i o be due 
to sulfur dioxide. Figure 2.6 illustrates the S0 2 cloud on April 6, approximately two days after 
the major eruption. This false color map of total ozone shows a region of brown colors surround- 
ing a black area extending east and west from the volcano located in southern Mexico. The black 
color indicates invalid ozone retrievals due to the S0 2 in the cloud. On this day the cloud was 
almost completely sheared into an eastward-drifting tropospheric component and a westward- 
drifting stratospheric component. The tiopospheric cloud was found to disappear within six 
days while the stratospheric cloud persisted and evolved slowly in shape due to wind shear. 

The El Chichon cloud was observed by three sensors on the SME spacecraft. The visible 
light spectrometer which observed the atmospheric limb radiance in the 400 nm region found 
enhanced scattering at latitudes near that of the volcano. Radiance perturbations were also 
found with the Near Infrared Airglow Spectrometer at 1.27 and 1.87 /am and the Infrared 
Radiometer at 6.3 and 9.6 /am. These sightings show the cloud to expand to cover the zone 
from the equator to about 30°N and to spread in altitude in the first three months after the 
eruption. 

2.2. 1.3 Airborne sampling surveys 

Research aircraft from the Aerosol Climatic Effects ( ACE) and Airstream Programs were 
deployed soon after the El Chichon eruption. While the less energetic Mount St. Helens plume 
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Figiue 2.6. TOMS total ozone map for 1230 LT, April 6, 1982 illustrating the apparent ozone excess due to S0 2 
absorption in the El Chichon cloud. The color scale for total amount of ozone is giver along the bottom. 



was very successfully sampled with aircraft, the primary cloud from El Chichon was about six 
kilometers above the peak aircraft flight altitude NVvt ther ss, valuable data were collected 
in secondary cloud layers and in the fallout from the primary cloud. Seven ACE sampling 
flights were conducted in April. May. and July 1982, and pre-eruption data were obtained in 
March. The ACE aircraft generally were equipped with the following instrumentation: 

1 . A osol samples 

• Knollenger Partick Size Spectrometer 

• Ames Wire Impactor 

• Quartz Crystal Microbalancv Cascade Impactor (QCM) 

• Condensation Nuclei Counter tl’niversity ol Minnesota) 

• Filter Collector 

2. Gaseous composition 

• Cryogenic cell 

• Lvman-o water vapor absorption cell 

The Airstream aircraft are equipped to obtain ambient air samples for laboratory 
analysis of trace gases and to collect aerosols with an impactor. 

2. 2. 1.4 Airborne lidar 

Information about the height and latitudinal extent of FI Chichon cloud was obtained 
during flights of the langley Research Center (LaRC) airborne lidar in July and October 198 



The July data, illustrated in Fig. 2.7, show that the primary cloud at 26 km had not spread 
beyond 30°N. In the second mission the cloud boundaries were at 3S°N and 10°S. 

2.2. 1 .5. Surface observations 

The El Chichon cloud fortuitously passed over the NOAA station at Mauna Loa about five 
days after the eruptions. This station is equipped to measure atmospheric spectral transmission, 
diffuse skylight, total ozone, and aerosol vertical distribution. The sensors used for these 
observations are: 

• NOAA Environmental Research Laboratory (ERL) Sun Photometer 

• Eppley Normal Incidence Photometer 

• Eppley Pyrheliometer 

• Dobson Spectrophotometer 

• Ruby lidar 

In addition to these a Canadian Atmospheric Environment Service — Brewer instrument was brought 
in after the eruption to measure column SO2 ■ 

Additional lidar observations are available from Fvkuoka, Japan (33.6S°N, 130.3S°E), Langley, 
VA (37. l°N, 76.3°W), L’Aguila, Italy (42.37°N, 13.4°E), and Garmisch-Partenkirchen, West 
Germany (47. 5°N, 11.0°E). 
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Figure 2.7. Langley Research Center lidar measurements of the ratio of the lidar to molecular back-scatter 
intensity of fil Chichon aerosol profile from 35°N to 13°N in July 1982. The arrows 

indicate the height of the tropopause. 



2.2.2. Cloud altitudes 


Some information about cloud altitude is available from wind data. Stratospheric winds were 
approaching the peak of the easterly phase of the quasi-biennial oscillation (QBO). while mid- 
tropospheric winds were westerly. The null point between these wind regimes from Hilo radiosondes 
was near 21 km. Thus a crude altitude discrimination is obtained from the drift directions of the 
cloud volumes. 

The northern part of the stratospheric cloud, according to TOMS and AVHRR analyses, 
drifted over Hawaii where detailed height data on the aerosol component was obtained with 
lidar soundings beginning at 0600 GMT, April 10. The primary layer was found at 25-26 km 
with secondary layers at 22 and 23.7 km. A persistent doud top at 26 km was found during 
succeeding observations on April 1 1 and 14. The lidar returns decreased to pre-eruption levels 
on April 23 and strengthened again after May 5 as the cloud presumably completed its drift 
around the world. These and succeeding soundings exhibited a peak near 27 km and a broader 
height distribution than on the initial transit. Lidar data from Fukuoka, Japan (33.65°N) 
beginning in May and the Langley lidar survey flight in July confirm that the layer with the 
greatest backscatter ratio is near 26 km. 

2.2.3. Cloud composition 

The primary cloud layer at 26 km could not be dirr ,t!y sampled from aircraft due to its height. 
However, lower altitude layers and fallout from the upper layer were sampled during Ames ACE 
and DOE Airstream flights. Particles collected with the Particle Size Spectrometer and QCM cascade 
impactor over Baja California on April 19 consisted primarily of acid-coated silicates peaking in 
number at 2 pm diameter with a total mass loading of 1-4 pg/m 3 . Ash particles collected on sbdes 
were found in large clusters. Sulfate particles collected by the Ames Wire Impactor were about 0.8 pm 
in diameter with a density of 0.6 pg/'m 3 . The highest sulfate concentrations (~ 160 parts per billion 
by mass (ppbm) as compared with a normal value of 1-2 ppbm) ever collected in 10 years of data 



were found during an Airstream flight on April 1 5. The concentrations of condensation nuclei 
was found to be larger than background levels at the edges of the cloud. 

Gaseous components, obtained on several Ames flights by a cryogenic collector, were analysed 
for SC^.COS, and C0 2 - All S0 2 mixing ratios were low (50 X 10“ 12 ) in April, although values 
greater than 100 X l(T i2 were found above 21.8 km in a later flight ir luly. The COS distribution 
was unusually high with 350 X 10" 12 found at *8 km on May 5. Airstream gas samples, however, 
showed no increase in COS up to 20 km. 

Water vapor was measured by Ly a absorption and showed an increase to 5 ppm compared with 
the long-term tropical average of 4.4 ± 0.4 ppm at 19-20 km. 

Column S0 2 was measured using Brewer spectrophotometers (Evans, 1982) at Mauna Loa 
beginning a few weeks after the eruption and on Ames CV-990 flights in May. These measurements 
generally indicate elevated S0 2 amounts of 10-16 m atm cm at Mauna Loa and 5-8 m atm cm at 
latitudes from 20-33°N during the Ames flight. However, pre-eruption measurements from Jamaica 
also indicated 8 m atm cm S0 2 . Thus, the S0 2 contritution by El Chichon is questionable in this 
data set. Later flights of this instrument on Airstream flights in July also found S0 2 column 
amounts of 10 m atm cm at Houston and 20 m atm cm near the Equator. 

The sulfur dioxide content of the initial cloud was measured with the Nimbus-7 TOMS, since 
S0 2 constitutes an interfering gas, which can be discriminated from ozone by the spatial signature 
of the cloud. Peak column S0 2 amounts during the first days after each mqjor eruption were about 
60 m atm cm. The total S0 2 content, obtained by integrating across the cloud, was 3-4 X 10 6 tons 
for the April 4 eruption. This S0 2 amount is considerably smaller than estimates derived from the 
Evans’ data. 
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2.2.4. Cloud optical properties 

The optical thickness of the El Chichdn cloud over Mauna Loa was measured using sun photom- 
eters (DeLuisi et aL, 1982) as shown in Fig. 2.8. In the initial passage of the cloud between April 9 
and 20, dust optical depths from 0.6 to 0.7S were found at 425 nm. The optical depth briefly 
dipped to near background levels (<0.1) on April 22-23, then returned to levels near 0.3 after 
April 27. The dust then slowly diminished to optical depths near 0.25 by July. The optical 
depths vary with wavelength, depending upon the chemical composition and size distribution, 
which change with time. King used this characteristic to deduce particle size distributions as 
discussed in section 3.2 

A. J. Krueger 
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2.3. SUMMARY OF SULFUR OXIDATION PROCESSES 

2.3.1. Introduction 

Sulfur compounds enter the atmosphere from anthropogenic and biological activities, sea salts, 
and volcanic emissions. On average, volcanoes may contribute about one third of the total sulfur 
emitted to the atmosphere, principally in the form of sulfur dioxide (Cadle, 1980). Local and 
regional concentrations of sulfur may greatly exceed those normally found in the atmosphere during 
and following volcanic eruptions. The constituents emitted by volcanoes are sulfur dioxide (SO 2 ), 
hydrogen sulfide (H 2 S), carbonyl sulfide (COS), and carbon disulfide (CS 2 ). Of these, SO 2 is 
normally the most abundant and is found in the emissions of all volcanoes. Table 2.1 from Turco 
et al. (1982) gives the concentrations of sulfur compounds measured in the Mount St. Helens’ plume 
and there is no reason to suspect that these relative amounts are other than typical. Once emitted 
into the atmosphere, all further chemical processes affecting sulfur are such as to produce higher 
oxidation states ending finally in the production of sulfate. The general nature of these processes 
may be classified as homogeneous gas-phase chemistry or as heterogeneous processes involving 
liquid phase or surface reactions. In no case is the complete oxidation sequence well understood, 
but several steps are either known or believed to be important. These will be summarized below. 

2.3.2. Gas phase oxidation of sulfur 

Gas phase chemistry provides the principal mechanisms for the oxidation of reduced 
sulfur to SO 2 and may contribute to the further oxidation of SO 2 to sulfate. Table 2.2 lists 
reactions of known or potential importance to these processes. The reactions R7 and R8 are 
functions of temperature, which is expressed in degrees Kelvin. Graedel in the reference column 
refers to Graedel’s (1977) review paper that cites the original references for the various reaction 
rates. 

The most highly reduced sulfur compound emitted from volcanoes is hydrogen sulfide. It is 
usually emitted in much smaller quantities than is S0 2 and is relatively rapidly oxidized in the 
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Table 2.2 

Gas Phase Oxidation of Sulfur Compounds 
(The numbers in parentheses are exponents to the base 10) 



Reaction 

Rate (cm 3 sec -1 ) 

Reference 

Rl) 

h 2 s + oh -♦ hs+h 2 o 

4.2 (-12) 

WMO 

R2) 

l 

HS + 0 2 -*> SO + OH 

<1 (-13) 

Graedel 

1 

HS + O -*> SO + H 

1.6 (-10) 

Graedel 


so+o 2 -»• so 2 + o 

2.7 (-17) 

Graedel 

R5) 

so + no 2 -> so 2 + no 

8.3 (-12) 

Graedel 

R6) 

CS 2 + OH -*• Products 

<1.5 (-15) 

WMO 

R7) 

C^ + O -*■ CS + SO 

2.4(-ll)Xe _530/T 

WMO 

R8) 

COS + 0 -*■ co + so 

2.1 (-11 ) Xe" 220 °/ T 

WMO 

R9) 

so 2 + oh -» hso 3 

6.7 (-13) 

Graedel 

RIO) 

so 2 + ho 2 -> so 3 +oh 

9.0 (-16) 

Graedel 

Rll) 

S0 2 + CH 3 0 2 -*■ SOj + CHjO 

<3.3 (-15) 

Graedel 

R12) 

so 3 + h 2 o -*■ h 2 so 4 

1.0 (-12) 

Graedel 


atmosphere. The lifetime against oxidation by OH is about a day. Once the thyil radical (HS) is 
formed it is in turn rapidly oxidized to sulfur monoxide (SO), reactions R2 and R3 of Table 2.2. 
The final step in the oxidation of H 2 S to S0 2 is oxidation of SO by 0 2 or N0 2 and this is essen- 
tially instantaneous. 

Carbonyl sulfide and carbon disulfide are much longer lived species than is hydrogen sulfide, 
though the oxidation rates are not yet well known. Both of these species may persist for times in 
excess of a month. It is unlikely that any of these reduced sulfur compounds are emitted in suf- 
ficient quantities to contribute significantly to the production of S0 2 and sulfate in a volcanic 
plume. 
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Sulfur dioxide is the predominant sulfur species emitted by volcanoes. Production of sulfuric 
acid and sulfates may occur to some extent through gas phase chemistry. The initial steps are oxida- 
tion of S0 2 by the odd-hydrogen radicals OH and H0 2 , reactions R9 and RIO of Table 2.2. There 
are additional possible gas phase reactions involving production and reaction of HSOx radicals, 
with x equal to 3-6 (Davis, 1979), but these are speculative and are not listed in Table 2.2. It is 
likely that heterogeneous processes play the major role in oxidation of S0 2 - 

2.3.3. Heterogeneous oxidation of S0 2 

Table 2.3 shows some of the liquid phase reactions involved in conversion of S0 2 to sulfate 
in droplets. The equilibrium constants for the reversible reactions shown in the table are the ratios 
of reactants to products, all expressed in units of moles/liter. The Adamowicz (1979) reference lists 
various aspects of the liquid phase oxidation of S0 2 . 


Table 2.3 

Liquid Phase Oxidation of Sulfur 
(The numbers in parentheses are exponents to the base 10) 



Reaction 

Equilibrium Constant 

Reference 

Rl) 

(S0 2 ) g 

(S0 2 )g 

3.3 (-2) 

Adamowicz 

R2) 

(S0 2 )g 

«-*■ H+ + HSOJ 

1.7 (-2) moles r 1 

Adamowicz 

R3) 

HSOj 

<-► H+ + SO 3 

6 2 (- 8 ) moles C -1 

Adamowicz 

R4) 

HSOJ (tm) 

- soj 



R5) 

SOJ (tm) 

- SO= 



R 6 ) 

h 2 so 4 

— H+ + HSOJ 

0.4 moles IT 1 

Adamowicz 

R7) 

HSO 4 

H+ + S0 4 

1 .2 (- 2 ) moles C -1 

Adamowicz 

R 8 ) 

HSO 3 + O 3 

-*• HSOJ + 0 2 


Turco et al. 

R9) 

HSO 3 + h 2 o 2 

- hsojo + h 2 o 


Turco et al. 


(tm) above refers to trace metal catalysis. 

The subscripts g and C on (S0 2 ) refer to gas and liquid phase concentrations. 
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An outstanding problem in this area is the mechanism by which trace metals catalyze the oxida- 
tion of bisulfite (HSOp and sulfite (SO 3 ). These processes are indicated schematically by R4) and 
RS) in Table 2.3. Another difficulty in describing the liquid phase oxidation of S0 2 is the fact 
that strong oxidants such as and, to a lesser extent, 0 3 are observed to significantly enhance 
the rate of production of sulfate (reactions R 8 to RIO). SO 2 absorption, and the subsequent 
equilibria of sulfur species in liquid drops are pH dependent. The presence of ammonia or other 
bases in the atmosphere may thus increase the capacity of liquid droplets for sulfur uptake. 

Reaction of S0 2 on the surfaces of suspended particles is another heterogeneous process 
influencing the rate of sulfur oxidation. The process is usually described in terms of a measured 
reaction efficiency for various solid materials. The details of the oxidation mechanism are not 
known and the terms representing this piocess in model calculations are empirical. 

Because of the dependence of heterogeneous oxidation of sulfur on factors such as the presence 
of trace metal catalysts, strong oxidants, base constituents, end the variety of solid surfaces present, 
a detailed theoretical description of particular episodes of S0 2 oxidation in the eruption plumes of 
volcanoes is difficult. A study of SO, oxidation in the Mt. St. Helens’ plume has recently been 
carried out by Turco et al. (1982). however. 


P W. Stewart 
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3.1. THEORY OF RADIATION 

The qualitative effects of the El Chichdn stratospheric cloud of aerosols and gases on radiances 
measured at a satellite can be deduced from a solution to the equation of radiative transfer: 

1(0 ,M,0) = I(t b ,M)e' tbSec0 

+ \ *<t ,X,3,0) 

o 

exp (-t sec 9 ) sec 0 dt, (1) 

where I (0, X, 9 , 0) is the specific intensity, or spectral radiance, of the radiant energy at *he 
aperture of a satellite radiometer; the first independent variable within parenthesis represents dis- 
tance above the c \..iace of the earth. It is expressed in terms of optical thickness (t), which increases 
from 0 at the top of the atmosphere to t b at the base of the atrnos P . • " other independent 

variables are X for wavelength, and 9 and 0 give the zenith and azimuthal angles, respectively, of the 
direction that the radiation is propagating. 

The first term on the right-hand side of Eq. (1) represents the intensity at the surface, but 
attenuated along the slant path to a satellite. The atmospheric residue from the volcano increases 
this attenuation, as well as altering the intensity of the sunlight transmitted to the surface to be 
subsequently reflected. The second term on the right-hand-side of Eq. (1) contains the source 
function (J), which accounts for radiant energy added along the line-of-sight. The volcanic residue 
in the atmosphere can contribute to J by both scattering and thermal emission processes 

Further insight in the nature of the radiation sources is gained by giving the complete equation 


for the source function: 
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J (t, X, 0, *) - (t, X) [JL (t, X, 0, 0, 0 O , 0 O ) X 

2n v 

S 0 (X) e Hsec + j f (t,X,0,0 5,0') X 

O *0 

I (r. X, 0', 0') sin 0* d0' d0*J + (1 - W 0 ) B(T(t), X), (2, 

where cj 0 is the albedo of single scattering, or the ratio of the scattering and extinction coefficients; , 

( 1 - w 0 ) gives the probability of absorption, and also for emission when the atmosphere is in local 
thermodynamic equilibrium. The dependent variables for the scattering phase function (p/4a) are I 

arranged such that the first pair of polar angles give the direction in which the energy is scattered j 

J 

and the last pair of angles give the direction of the incident pencil of energy. B is the Planck fune- j 

i ' 

tion at level t where the air temperature is T. 

] 

The terms within the square brackets (times u> 0 ) of Eq. (2) account for the scattered energy. 

The first of these terms gives the contribution of the direct sunlight, and the other term gives the 
contribution of the diffuse radiant energy arriving from all directions (0', 0') and scattered in the 
direction(0, 0). The last term (outside of the brackets) accounts for the thermal emission. 

A month or so after the last El Chichon eruption, lidar scattering trom the aerosols indicated 
that the El Chichon cloud was confined to the layer between 1 8 and 35 km, with peak concentration 
near 26 km. The trace gases from the volcano are of such low concentration that their scattering 
effects can be neglected. The gases absorb and emit energy, however. As will be shown in section 
3.2, the aerosol particles are of submicron size. They strongly scatter light in the visible and ultra* 
violet spectrum; but his effect decreases in the infrared, while their absorption increases. 

The vertical profile of the source function depends on the profiles of the absorbing and scatter- 
ing constituents. If under normal conditions the source function peaks high in the stratosphere 
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and is zero below 35 km, then the volcanic cloud does not affect the radiant energy reaching a 
satellite. Such is the case for X < 290 nm in the SBUV experiment. Generally, however, the source 
functions indicate a stratospheric response at or well below 26 km under normal conditions (see 
Figs. 4.4 to 4.6). Then the volcanic layer will modify the source function and affect the radiant 
energy coming from the earth and atmosphere. The magnitude of the effect depends on wavelength 
and will be discussed in the following sections. 


R. S. Fraser 

3.2. RADIATIVE CHARACTERISTICS OF THE AEROSOLS 

The effect of aerosols from the El Chichon eruptions on satellite remote sensing depends 
primarily on their optical thickness, phase function, and refractive index. After a few weeks these 
aerosols are restricted to the stratosphere. For some remote sensing experiments, sveh as limb 
scanners or solar occultation observations the vertical distribution of the aerosol layer will also play 
a key role in modifying the radiation field at satellite altitude. 

In order to determine the aerosol optical thickness and aerosol size distribution, a series of 
solar radiation measurements from Mauna Loa Observatory was collected (DeLuisi et al., 1982). 
The monthly mean stratospheric aerosol optical depths were determined from these data for June 
and July by subtracting the molecular scattering, ozone absorption and background tropospheric 
aerosol contributions from the total optical depth. These results, presented in the left portion of 
Fig. 3.1, show that the stratospheric aerosol layer resulting from the El Chichon eruptions had an 
optical depth between 0.15 and 0.25 throughout the visible wavelength region. This is more than 
an order of magnitude larger than background levels at Mauna Loa (Shaw, 1979) and larger than 
all but a small percentage of days at Tucson (King et al.. 1980). 
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The columnar aerosol size distribution can be determined from the spectral variation of the 
aerosol optical depth using constrained linear inversion methods (King et al. 1978; King. 1982). 

The size distributions obtained by inverting the Mauna Loa measurements are presented in the right 
portion of Fig. 3. 1 . while the solid curves in the left portion indicate how the inverted size distribu- 
tions are able to reproduce the aerosol optical thickness measurements. Due to the prominent 
negative curvature in the spectral optical thickness measurements, the corresponding size dtetribu- 
tions are quite narrow with peak concentrations at radii between 0.3 and 0.4 pm. Though all 
inversions were performed assuming the complex refractive index of the aerosol particles was wave- 
length and size independent and given by m * 1.45-0.00i, it is well known that the refractive index 
sensitivity is quite weak, affecting the inverted size distribution slightly by shifting its magnitude 
and radii while maintaining its overall shape (King et al.. 1978) 

The columnar size distributions presented in Fig. 3. 1 repres the number of particles per 
square cm per log radius interval in a vertical column through the stratosphere. For a particle density 
of 2 g enf \ the columnar mass loading for these cases is 0.0S79 g m“* and 0.0602 g m"* for June 
and July, respectively. These mass loadings were obtained over the radii range 0.1 < r < 4.0 pm 
by extrapolating the inverted size distributions to smaller and larger radii. 

Both the June and July size distributions are well described by a modified gamma distribution 
of the form 

= Cr° exp(-br). (1) 

dlogr 

where the mode radius r 0 = a/b. Each distribution presented in Fig. 3. 1 was fit to Eq. ( 1) using 
an efficient gradient-expansion method from nonlinear least-squares theory (Bevington, 1969). 

For the June distribution the best-fit parameters were C = 1 .549 X 10**. a s 8.26. b = 24.4 pm“’ , 
r 0 = 0.338 pm and for July C = 9.897 X I0 ,g . a = 1 3.65. b = 39.3 pm' 1 . r Q = 0.347 pm. 
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In addition to the mass loading, the inverted size distributions can be used to determine the 
optical thickness and phase matrix elements as a function of scattering angle and wavelength. 
Assuming the aerosol particles to be composed of H 2 S0 4 * bH 2 C droplets having a 75% concentra- 
tion by weight of sulfuric add, we obtained the aerosol optical thickness, single scattering albedo 
and asymmetry factor for wavelengths between 0.25 and 25.0 pm. For these computations, 
presented in Figs. 3, 2-3. 4, we used the July size distribution presented in Fig. 3.1 together with the 
complex refractive indices tabulated by Palmer and Williams (1975). Similar computations were 
performed using the modified gamma distribution having the appropriate parameters for the July 
distribution, but the results were virtually identical to those presented in Figs. 3.2-3.4. The optical 
thickness and single scattering albedo are the most important parameters governing the radiative 
properties of the aerosol layer in the near and thermal infrared. With the values of the optical thick- 
ness and single scattering albedo from Figs. 3.2 and 3.3 for X = 3.7, 10.8 and 12.0 pm, for example, 
the El Chichon layer would be expected to produce a -0.6K bias in the AVHRR derived sea surface 
temperature at night (-1 .2K during the day). In the 1 1.0 pm window region the emissivity of the 
layer is nearly unity (e = 1 - w 0 ) but the optical thickness has a strong wavelength dependence in 
this region. At 3.7 pm the emissivity is less and the aerosol optical thickness is greater than at 1 1 pm. 

To examine the phase matrix in further detail, we computed the phase function and degree of 
polarization of the aerosol particles composing the stratospheric layer for a wavelength X = 0.7 pm. 
This phase function, illustrated in Fig. 3.5, has an asymmetry factor g = 0.768, a single scattering 
albedo <o 0 = 1 .000, and a ratio of the volume extinction to volume back-scattering coefficients 
S = 73.3 sr. The latter parameter has been found to be especially useful in analyzing monostatic 
lidar observations (Spinhime et al., 1980; Reagan et al., 1980). The degree of polarization, illus- 
trated in Fig. 3.6, has an unusually large number of scattering angles at which the polarization is 
neutral. 
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Figure 3.2. Aerosol optical thickness as a function of wavelength for the July 1982 
stratospheric aerosol size distribution presented in Figure 3.1, where the aerosol 
particles are assumed to be composed of a 75% H 2 S0 4 solution. 
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Figure 3.4. As in Figure 3.2 except for the asymmetry factor 
of the aerosol particles. 
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Figure 3.5. Phase function as a function of scattering angle for the July 1982 
stratospheric aerosol size distribution presented in Figure 3.1, where X = 0.7 jum 
and m = 1.428 - 2.07 X 1 0~®i. This phase function has an asymmetry factor 
g = 0.768 and an extinction to backscatter ratio S = 73.3 sr. 
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Figure 3.6. Degree of polarization (I r - Ig)/(I r + Ig) as a function of 
scattering angle for incident unpolarized radiation and for the July 1982 
stratospheric aerosol size distribution presented in Figure 3.1, where 
X = 0.7 pm and m = 1.428 - 2.07 X 10" 8 i. 
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During the forthcoming months radiative transfer computations will be performed of the 
transmitted intensity and degree of polarization in the principal plane of the sun. These computa- 
tions will initially be made using the optical thickness and phase matrix results presented in Figs. 3.1, 
3.5 and 3.6 for July 1982 Mauna Loa observations. In addition, the vertical distribution of the 
aerosol particles obtained from monostatic lidar will be used. These computations will be compared 
with surface measurements collected b. Coulson (1982) to validate our radiative transfer model and 
to verify that we understand the optical properties of the aerosol layer. Once satisfactory agreement 
has been achieved, the radiative transfer model will be used to examine the expected effects of the 
stratospheric aerosol layer on satellite observations. This research will be a collaborative effort 
between R. S. Fraser and M. D. King. 

In addition to this investigation, a series of spectral optical thickness measurements between 
0.44 and 0.87 jum have been collected by Dr. J. D. Spinhime during a three week period in October 
and early November 1982. The solar transmission measurements were collected from NASA’s 
Electra aircraft which flew from Wallops Island, VA to Santiago, Chile, and beyond. This flight 
covered the latitude range from 49°N to 46° S and data from the cloud-free daylight portion of the 
flight will be used to derive the latitudinal distribution of total ozone, aerosol optical thickness and 
aerosol size distribution. 


M. D. King 
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4. EFFECTS OF EL CHICHON ERUPTIONS ON SATELLITE OBSERVATIONS 


4.1. NOAA/AVHRR PRODUCTS 

4.1.1. Change in albedo 

Satellite measurements of planetary radiation Lorn NOAA-7 AVHRR data have been scruti- 
nized to examine the impact of El Chichon. It was known in May 1982 that an overall increase in 
albedo was produced because the number of daytime satellite retrievals of sea surface temperature 
(SST) (see below) had decreased markedly in the latitude band affected by El Chichon. This 
decrease has been shown to be due almost entirely to a test using a threshold for reflected radiation. 
This test does not accept any SST retrievals unless the albedo in AVHRR channel 2 (X 0.71-0.99 pm) 
falls below 1 .7%. This albedo is referenced to the response that would occm by the AVHRR observ- 
ing a Lambert surface with 100% albedo, when exposed to the solar radiant energy above the atmos- 
phere. 

Albedo data from the NESDIS operational pk netary radiation budget archive have been 
examined using a 10-day minimum albedo classification. Minimum daily albedos at every 2.5° 
lat/long grid were retained globally over each 10-day period. A composite minimum albedo was 
constructed every month (3-ten day periods) using a simple arithmetic mean of all three values at 
each grid location. These may be contrasted with a similar sequence in 1979. The albedos for the 
areas of minimum cloudiness between 0° to 30°N are greater in 1982 than in 1979. This anomaly 
was still persisting during September 1982. 

Four target areas within this aerosol belt were selected for further examination. Measurement 
of albedo in cloud-free areas was noted immediately prior to the passage of the dust cloud (see 
section 2.1) and just after dust cloud passage. In this table the albedo was corrected for the cosine 
of the solar zenith angle. Cloud-free scenes were examined further for daytime outgoing longwave 
radiation. The results are shown in Table 4.1. Albedos over the cloud-free ocean have nearly 
doubled, from about 5% to 10%, while outgoing flux shows a several watt/nri decrease. 


4-1 



SSSSffl* 


Table 4.1 

Cloud-Free Radiation Changes for Four Areas: Eastern Pacific (10-20N, 10" 1 !SW); 
Philippine Sea (10-20N, 125-135E); Bay of Bengal (5-15N. 85-95E); 
and Arabian Sea (10-20N, 60-70E). 



Albedo (%) 

I-icrease (%) 

Outgoing Long Wave 
(11-12 pm; WM- 2 ) 

Decrease 

'WM-2) 

Eastern Pacific 





2 April 

5.1 


304 


6 April 

10.4 

5.3 

294 

10 

Philippine Sea 





14 April 

6.0 


299 


15 April 

7.3 

1.3 

296 

3 

Bay of Bengal 





6 April 

4.5 


305 


18 April 

10.5 

6.0 

291 

14 

Arabian Sea 





1 5 April 



303 

1 

19 April 

u 

6.3 

296 

i_J 


4.1.2. Sea surface temperature 

As mentioned above, AVHRR processing for multichannel sea surface temperature (MCSST) 
initially responded to the El Chichon eruptions with a drastic reduction in daytime retrievals and a 
pronounced negative bias (temperature lower than those from ships) on those nighttime observations 
that passed all cloud tests. Drifting buoy validation statistics had been showing ±0.5°C accuracies 
(RMSD^ before the eruption and continue to show similar [±0.5°C] accuracies outside the contam- 
inated belt. Initial biasing was seen at fixed buoy matchup positions within the aerosol cloud at 
locations such as Hawaii and in the Gulf of Mexico. One such exam le is given in Fig. 4.1 using 
differences from a large moored NOAA buoy (#42003) in the central Gulf of Mexico. The final 
eruption on 4 April (Julian Day 94) marks the start of a -1.5°C bias that remained virtually constant 
until 25 J *. when apparently another level of aerosol moved north over the buoy location causing 
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Gulf of Mexico (26N, 86W) during a portion of 1982. 


a greater than -2°C offset in derived MCSST from AVHRR. The absolute value of the bias decreased 
after the summer by 2°C, perhaps as a result of ihe stratospheric aerosol layer thinning out. 

The final series of Fig. 4.2a-e shows the apparent extent of the volcanic cloud on a monthly 
mean basis using this negative bias relative to ships as an indicator of aerosol intensity. Monthly 
mean data were differenced using ship and buoy only and satellite only. Areas exhibiting differences 
exceeding 1C have been highlighted thereby revealing spread and diffusion of the cloud throughout 
the summer. 

The following observations can be made: 

1. Maximum contamination occurred during July, probably when SO 2 conversion to H 2 SO 4 
was nearly complete and before the larger particles had fallen out of the stratosphere. 

2. Some transport north of 30N has been observed on a monthly basis, especially over the 
central Pacific, but most of the aerosol continues to occvpy the latitudinal belt of 0 to 30N. 

3. Evidence of movement into the Southern Hemisphere t o 1 OS has been observed since 
June, particularly over the South Equatorial Atlantic Ocean. (The paucity of ship observa- 
tions may be responsible for holding back any spread of negative bias in the MCSSTs in the 
central and eastern South Equatorial Pacific.) 

4. Beginning in August the area around Indonesia became incorporated into the El Chichon 
belt. It is possible that the June-August eruptions of Galunggung were responsible for this. 

5. During September an apparent migration of aerosol northward over Japan has been noted. 
Since August, lidar observations from Fukuoka, Japan have been showing more aerosol, 
particularly at 20-23 km (below El Chichon ’s 25-30 km altitude), that may be moving 
northward from Galunggung in Indonesia. 
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Figure 4.2. (Continued) 



6. By September the general levels of bias exhibited by the satellite-derived MCSSTs appear 
to have lessened over the Western Hemisphere, although the distribution has remained 
largely unchanged since July and August. Lower levels are presumed to be attributed to 
fallout of larger particles from the stratosphere, which emit radiation at lower temperatures. 

The negative bias in the MCSST sea surface temperature decreases after the 230 day (Fig. 4.1). 
Perhaps this change results from smaller (<1 micron) particles remaining after the larger ones 
(>1 micron) fall out into the troposphere. We are investigating several aerosol correction schemes at 
NESDIS for incorporation into the MCSST procedures. If successful, this would provide a means for 
continuous monitoring of aerosol, as well as providing a mechanism to remove subsequent volcanic 
contamination effects on the retrieval of sea surface temperatures from satellite inf.a-ed measure- 
ments. 

4.1.3. Aerosol remote sensing 

An algorithm has been developed and tested in ground-truth experiments to convert the 
radiance measured by AVHRR Channel 1 (0.56-0.68 nm) over cloudless ocean areas into an aerosol 
optical thickness at a wavelength of 0.5 /urn (Griggs, 1 98 1 ). A table look-up procedure is used. The table 
was constructed using the multiple scattering code of Dave for a plane parallel atmosphere with 
Lambertian Earth surface reflectance. The parameters of the model were chosen to give agreement 
between theory and observation using Landsat 1 MSS data and Volz type sun-photometer solar 
extinction measurements off the coast of San Diego. Figure 4.3a shows a comparison of the optical 
thickness derived from this algorithm applied to NOAA-6 AVHRR Ch 1 data and sun-photometer 
extinction measurements at several locations used in a 1980 groundtruth experiment. One N unit 
is equivalent to an aerosol optical thickness of 0.21 3 at X 0.5 /am. The two estimates are in good 
agreement. The correlation coefficient is 0.93 and the standard error about a linear regression line 
is ±0.25N (±0.05 in optical thickness). 
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Figure 4.3b shows a comparison between aerosol optical thickness derived from NOAA-7 
AVHRR data and from sun-photometer data at the Mauna Loa observatory between March 29, 
1982 and April 28, 1982. The aerosol was first observed at Mauna Loa on April 9th. The area 
being studied with AVHRR data is about 800 km WNW of Mauna Loa. so the aerosol probably 
arrived in this area on the 11 th. However, due to constraints imposed by clouds or sun-glint 
effects, the first possible WHRR aerosol measurement occurred on the 16th. Since this area 
is quite distant from Mauna Loa and the AVHRR measures total optical depth to the sea 
surface, whereas the Mauna observations are made 3.4 km above the sea, the observations are 
in very good agreement. Of particular interest is the decrease in optical depth after the 18th 
at both locations, suggesting that the aerosol cloud had both a leading and trailing edge. The 
numbers next to the AVHRR poin>* represent the number of water cloud screening tests passed 
out of a total of eight. These tests are used in the MCSST program and indicate that the aerosol 
cloud is being classified as a water cloud by many of these tests; hence, no surface temperatures 
would be derived on these dates. 

A method is currently underdevelopment which uses Ch 1 and Ch 2 AVHRR data. 
Preliminary results suggest that information on aerosol size distribution as well as optical 
thickness can be extracted from a two channel algorithm. 


A. E. Strong 
A. Gruber 
L. Stowe 
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Figure 4.3b. El Chichon aerosol optical thickness versus time. 
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4.2 HIRS/MSU SOUNDERS 
4.2.1. Background 


HIRS and MSU are components of the operational infrared and microwave temperature and 
humidity sounders on the satellites NOAA-6 and NOAA-7. Table 4.2 shows the central frequency 
(i>) for each channel; the third column gives the pressure P at the peak of the weighting function, 
dr/dlnp, for each of the temperature sounding channels; r is the vertical transmission of the 
atmosphere. TV pressure corresponds to the level of the atmosphere where the typical atmospheric 
optical depth is one above that level. The fourth column gives the pressure level at the peak of the 
product of the weighting function and the Planck black body function (B) of the atmospheric 
temperature (T). This product, after being normalized to the total radiance (I) at the top atmos- 
phere, is plotted along the abscissa as a function of pressure for the HIRS and MSU channels in 
Figs. 4.4-4.6. The HIRS observations have a horizontal resolution of 18 km at nadir, changing to 
60 km at a 55° angle from the nadir. Under the same conditions, the MSU resolution changes from 
100 km to 300 km. 

Global atmospheric temperature retrievals are produced operationally from the HIRS2/MSU 
data. In addition. NASA has developed a research oriented processing system for the HIRS2/MSU 
data (Susskind et al., 1982) which deduces not only atmospheric temperature profiles but also sea 
surface temperature, ice and snow cover, and percent cloud cover and cloud height from this data. 
Two distinct questions are of interest. The first question is the extent, if any, that the quality of 
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Table 4.2 

HIRS2 and MSU Channels 


Channel 

v fern -1 ) 

Peak of 
dr/dlnp (mb) 

Peak of 
Bdr/dlnp (mb) 

HI 

668.40 

30 

20 

H2 

679.20 

60 

50 

H3 

691.10 

100 

100 

H4 

703.60 

280 

360 

H5 

716.10 

475 

575 

H6 

732.40 

725 

875 

H7 

748.30 

Surface 

Surface 

H8 

897.70 

Window, sensitive to water vapor 

H9 

1027.90 

Window, sensitive to O3 

H10 

1217.10 

Lower tropospheric water vapor i 

HI 1 

1363.70 

Middle tropospheric water vapor 

HI 2 

1484.40 

Upper tropospheric water vapor 

H13 

2190.40 

Surface 

Surface 

H14 

2212.60 

650 

Surface 

H15 

2240.10 

340 

675 

H16 

2276.30 

170 

425 

H17 

2310.70 

15 

2 

H18 

2512.00 

Window, sensitive to solar radiation 

H19 

2671.80 

Window, sensitive to solar radiation 

Ml 

50.30a 

Window, sensitive to surface emissivity 

M2 

53.74a 

500 


M3 

54.96a 

300 


M4 

57.95a 

70 



a values in GHz 


the retrieved products is degraded as a re«"'t of the volcano. The second question is the extent to 
.vhich sounding data can give us additional information about the volcanic eruption and its radiative 
effects. 


Thus far, the oyproducts of El Chichon eruption have been observed indirectly, as degraded 
products of ozone and sea surface temperature from analyses of TOMS and AVHRR data, and 
directly, from the visible imagery of AVHRR. It is thought that TOMS is observing primarily the 
effects of S0 2 gas from the eruption and AVHRR is being affected by aerosols. 
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Figure 4.4. Normalized source function for the MS J temperature sounding channels. 
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CHANNEL 

FREQUENCIES, cm 

HI 

668.60 

H2 

679.05 

rl3 

689.70 

H4 

703.80 

H5 

716.70 

H6 

731.85 

H 7 

749.50 

U.S. STANDARD ATMOSPHERE 
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Figure 4.5. Normalized source function for the H1RS 15 jum 
temperature sounding channels. 
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Figure 4.6. Normalized source function for the HIRS 4.3 /urn 
temperature sounding channels. 





While SO 2 affects observations in the channel of TOMS, SO 2 has no appreciable absorption 
in any of the channels of AVHRR or the temperature sounding channels of HIRS. Therefore, one 
would not expect it to degrade the retrieved temperature fields in any way. S0 2 does absorb, to 
some extent, in HIRS channels 11,10 and 18 in decreasing order of strength. The first two channels 
are primari'y for humidity sounding. The retrieved humidities may therefore be affected. The last 
channel, 18, is a shortwave window channel used by GLAS to obtain sea surface temperature. 

During the day, channels 18 and 19 are used together to solve simultaneously for the surface tem- 
perature and the bi-directional reflectance of solar radiation. At night, independent estimates of 
surface temperatures can be obtained from the radiances of each channel ( 18 or 19). In general, 
these estimates differ by less than 1°C. The difference in sea surface temperature retrieved from 
these two channels can serve as a tracer for S0 2 , The effect is expected to be very small however. 

The H 2 S0 4 aerosol is expected to have b»oad banded attenuation features throughout the 
3-1 5 jum regions sampled by the HIRS and AVHRR instruments. This attenuation is due almost 
exclusively to absorption, as shown in Fig. 3.2 and 3.3. Channels 3 (3.55-3.93 jum) and 4 (10.5-1 1.5 
jum) on AVHRR both contain frequencies characterized by relatively large H 2 S0 4 absorption and 
the effect of this absorption has been noted on the sea surface temperature fields retrieved from 
AVHRR. utilizing these channels. Sea surface temperatures from HIRS2 are retrieved by GLAS 
primarily from the two shortwave window channels 1 8 (3.95-4.00 jum) and 1 9 (3.7 1 -3.83 jum) which 
have considerably less absorption due to H-,S0 4 aerosol than AVHRR channel 3. The longwave 
window channel 8 (10.9-1 1.3 jum) has comparable absorption to that of AVHRR channel 4 but this 
channel is not used by GLAS to retrieve sea-surface temperature. 


It can be inferred from Komfield and Susskind (1977) that the brightness temperature in a 
window channel at frequency v(cm _1 ) will differ from the surface temperature Tg by approximately 

T|d-r> 


AT = 


1.439v 


( 1 ) 
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if the radiation is weakly attenuated. This implies that absorption at 1 1 pm is three times more 
effective than the same absorption at 3.7 pm in depressing the brightness temperatures. From 
Fig. 3.2, one would expect a brightness temperature error of about 0.8°C at 1 1 pm at nadir and 
3O0°K surface temperature, and likewise about 0.8°C in AVHRR channel 3 at about 3.7 pm. This 
error is magnified in the analysis of the operational AVHRR data because of the use of regression. 

In HIRS channels 18 and 19, however, the error is expected to be much smaller and potentially 
insignificant. In addition, in the physical retrieval method, no regression is used and an error in 
computed brightness temperature results in the same error in sea-surface temperature. 

An obvious after effect of the volcanic eruption is the visible cloud tracked around the globe 
by Matson (see section 2.1). This same cloud is readily apparent in most IR channels. By observing 
this cloud in a number of spectral channels sounding different portions of the atmosphere in dif- 
ferent spectral intervals, one can infer some of the height and optical thickness properties of the 
cloud. The following section shows preliminary results based on observations of the cloud from 
HIRS2 data on March 29. 10 GMT, 5 hours after the initial eruption, and March 30. 9 GMT. 

4.2.2. Preliminary observations 

From imagery of the infrared brightness temperatures for March 29 ’ 9 GMT. features associ- 
ated with the eruption are observed in all the HIRS2 channels but channe* . md 1 7, which sound 
the mid-upper stratosphere. In channel 8. the 1 1 pm window, an extremely cold cloud is observed. 
The coldest spot is at 17.07°N, 93. 14°W about 2.3° East of El Chichon, with a brightness 
temperature of 195. 1°K. The cloud runs east-northeast, west-southwest with its maximum north- 
east extent given at 18.35°N, 91.59°W, with a 229.6° brightness temperature, in channel 8, and 
maximum southwest extent given at 17.04°N, 95.40°W, almost coincident with El Chichon. with a 
239.6° brightness temperature. Outside the cloud, the biightness temperature rises sharply to a 
typical value of 286.7°, which was taken at 1 7. 19°N, 96.46°W. The brightness temperatures can be 
compared with the air temperatures measured at Ver ,iz. Mexico.at 12 GMT, 29 March (Fig. 4.7). 
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Brightness temperatures for all channels sensitive to the volcanic cloud are shown in Table 4.3, 
for both March 29 and March 30. On March 29, all observations are shown for the coldest spot (a), 
the southwest comer of the cloud located over the volcano (b), and a typical spot (c). In addition, 
observations in two window channels are shown for the northeast comer of the cloud (d). 

On March 30, the cbseivations are shown for equivalent locations marked e-h. It is significant to 
note that on March 29 signals due to the cloud are clearly observed in stratospheric sounding chan- 
nels 2 and 3, with the southwest edge of the cloud, located over the volcano, appearing warmer than 
the remainder of the field. Channel 2 is totally insensitive to observations lower in the atmosphere 
t!ian 200 mb (12 km), indicating that the cloud must be at least in the upper troposphere. 

If a cloud covers a spot with fractional cloud cover a, the radiance in channel i is given by 

I 4 = (1 -a) Ij,CLR + aI i,CLD 

where Ij CLR is the radiance which would be seen if the field of view were clear and Ij CLD is the 
radiance which would be observed if the cloud completely covered the field of view of the instru- 
ment. The cloud blocks contributions to the radiance emitted from beneath it which would other- 
wise be observed by the channel. In general, clouds are in the troposphere and appear colder than 
the air beneath them. Therefore clouds in general appear cold. The warm spot over the volcano 
seen in channels 2 and 3 implies that if the cloud is at the same temperature as the atmosphere, 
either the cloud immediately above the volcano is higher in the stratosphere than the surrounding 
clouds, which would increase Ij cld> or the a * r * n the stratosphere above the volcano is warmer 
than its ambient values, which would increase both Ij anc * Ij CLD- 

The coldest spot appears increasingly colder as we go from channel 2 to 8, in the order of 
decreasing atmospheric optical depth from the satellite to the cloud. This is consistent with mixing 
' i decreasing amounts of the warm stratosphere above the cloud with increasing spectral frequency. 
The brightness temperatures in the 4.3-3. 7 pm region, seen in channels 13-19, are consistently 
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Table 4.3 

HIIIS Brightness Temperatures 


CH 

March 29 
Coldest 8 

SW 1 * 

Typical® 

NE d 

March 30 
Coldest® 

SW f 

Typical* 

NE h 

1 



232.3 






E 

217.0 

221.6 

218.1 


217.4 


218.0 


E 

216.9 

221.3 

217.5 

i 

214.2 

216.3 

217.4 


H 

202.3 

218.9 

231.3 


216.8 

223.5 

231.0 


fl 

201.3 

228.6 

245.0 


220.7 


246.7 


B 

198.5 

231.6 

254.1 


223.6 


259.4 


7 

196.6 

235.4 

263.7 


228.4 


272.4 


8 

195.1 

239.6 

286.7 

229.6 

242.4 

248.8 

287.0 

268.0 

9 

213 6 

228.8 

266.9 


247.3 


273.1 


10 

213.3 

250.2 

279.7 


252.6 


286.6 


11 

197.3 

231.9 

251.7 


232.3 


261.6 


12 

197.8 

236.5 

248.8 


222.4 


241.4 


13 

205.1 

255.0 

275.6 


258.0 


283.2 


14 

205.6 

246.4 

266.3 


248.2 


269.9 


15 

210.1 

234.1 

250.9 


233.9 


251.2 


16 

17 

18 

217.6 

224.3 

230.3 


222.5 


228.8 


208.6 

262.0 

282.6 

232.9 

273.9 

281.5 

291.1 

282.6 

19 

212.1 

265.6 

283.0 


277.6 

282.3 

290.9 



a) 17.07°N 

93. 1 4° W 

e) 23.37°N 

89.42°W 

b) 17.04°N 

95.40°W 

0 22.28°N 

92.41°W 

c) 17.19°N 

96.46°W 

g) 21.41°N 

94.50°W 

d) 18.35°N 

91.59°W 

h) 23.94°N 

88.17°W 
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warmer than in the IS /im region. This implies that the cloud has bulk emissivity (e) less than 1 
at least in some spectral region. If ae were equal to 1, both the long and short wave window chan- 
nels should have brightness temperatures nearly equal to the temperature of the cloud. The sum of 
e plus p, the reflectivity of the cloud, +t, the transmissivity of the cloud equals 1. For non-unit e, 
the cloud radiance consists of e of the black cloud radiance, plus p of the very cold reflected 
space radiance, plus r of the very warm transmitted clear column radiance. If e is less than 1, then 
the brightness temperatures associated with the shortwave channels would always appear warmer 
than equivalent longwave temperatures because of mixing of the apparent brightness temperature 
from mixed radiance fields is weighted more toward the warm brightness temperature at the higher 
frequencies. This is consistent with the higher brightness temperatures at high frequencies. In ad- 
dition, the higher brightness temperatures at high frequencies may imply that the cloud has higher 
emissivity in the long wave region than in the short wave region. The same phenomenon is observed 
in the spot above the volcano to even a larger degree. 

On March 30, about 24 hours later, the cloud has approximately the same shape but has moved 
about 5°N and 3°E from its location on the day of the eruption. El Chichon appears to be cloud 
free at this time. The coldest spot has likewise translated by about the same amount. The appear- 
ance of the cloud in the infrared channels is quite different from the previous day, however. No 
warm spot exists in channels 2 and 3 and the cloud can barely be observed in channel 2. On the 
other hand, the signal of the coldest spot of the cloud is larger in channel 3 than in the previous day. 

The cloud effects in the other longwave cl anneis, 4-8, are clearly much less than in the previous 
day. Not only is the effect of the cloud smaller in each channel than on the previous day, but the 
brightness temperature in the coldest spot increases with decreasing optical depth, which is counter 
to the trend of the previous day. Because of the extreme cold temperature in channel 8 on March 
29, it is almost certain that the coldest spot was completely covered by an optically thick cloud 
for the long wave channels. On the other hand, on March 30, a(l -t) must be considerably less 
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than 1 for the long wave channels, even in the coldest spot. The short wave channels 13-16 and 18, 
19 show the same trend. In addition, the brightness temperatures in the short wave channels in the 
cloudiest spot are much warmer than in equivalent long wave channels. This is consistent with 
at(l - t) considerably less than 1, and in addition, with r, the transmissivity, increasing with decreas- 
ing wavelength. 

No signal attributable to El Chichon was observed in any MSU channel on either day. This may 
be in part due to the lower spatial resolution of the channels. In addition, no apparent signal was 
observed in any HIRS2 channel in the next satellite pass, March 30, 20 GMT. 

4.2.3. Future plans 

More comprehensive studies of the effects of El Chichon ejecta on retrievals of atmospheric 
temperature, sea surface temperature, and cloud fields will be done for a long period of time com- 
mencing with March 1982 in order to assess the effects of volcanic eruption on retrieved atmospheric 
temperature profiles and sea surface temperatures and to track the volcanic cloud and obtain its 
properties. Of primary interest is whether accurate sea-surface temperatures can be obtained from 
HIRS2, even under the H 2 S0 4 aerosol cloud caused by the volcano. 

J. Susskind 


References (Section 4.2) 

Komfield, J., and J. Susskind, 1977: On the effect of surface emissivity on temperature retrievals. 
Mon. Wea. Rev., 105, 1605-1608. 

Susskind, J., J. Rosenfield, D. Reuter, and M. T. Chahine, 1982: The GLAS Physic-*! Inversion 
Method for Analysis of H1RS2/MSU Data. NASA Technical Memorandum 84936. Goddard 
Space Flight Center, Greenbelt, MD 20771. 


4-24 



4.3. VASONGOES 


4.3.1. Expected effects 

The VISSR Atmospheric Sounder (VAS) is a 1 2-channel radiometer located on the operational 
geosynchronous GOES satellites. Twelve VAS channels are located in the 15 and 4.3 pm C0 2 bands, 
the 6.7 and 12 pm H 2 0 bands, and the 1 1 and 3.9 pm windows. Consequently, the spectral effects 
of El Chichon on the VAS channels will be the same as on the corresponding HIRS channels (Sec- 
tion 4.2). VAS data is archived routinely at the University of Wisconsin (UW) and occasionally by 
the VAS Demonstration at GSFC. NOAA’s operational VISSR schedule rarely permits 12-channel 
VAS data collection from low latitudes. VISSR operations do permit hourly 3-channel scans of the 
globe in the 1 1 and 12 pm split window channels, coupled with either the 3.9 pm window or the 
6.7 pm upper-air water vanor channel. 

a) The VAS split window channels are used to estimate sea surface temperature (SST). 

An empirical regression formula developed by UW for SST along the East Coast of the 
United States is: 

SST » 2.23 Tj* - 1.23 T |* 2 , (1) 

where Tj* and T * 2 are the apparent temperatures in the 1 1 pm and 12 pm channels, 
respectively. The effects of El Chichon on VAS SST should be detectable at the ±1°C 
level, provided there are no other clouds in VAS’ 15 km footprint. At low latitudes, the 
3.9 pm window is much too affected by reflected sunlight to be useful during daylight 
operations. 

b) The VAS 6.7 pm water vapor channel is like the familiar METEOSAT band, with its 
radiance representing the average brightness of the top few millimeters of the atmospheric 
water vapor. The 6.7 pm radiance normally arises from a layer a f, w hundred millibars 
thick, centered at a level between 300 and 500 mb. However, molecular water vapor 
injected by El Chichon into the lower stratosphere would appear as an unusually cold 
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feature in the 6.7 pm radiance field. For example, 1 mm of predpitable water in the layer 
from SO to 250 mb (mixing ratio about 0. 1 gm/kg, or relative humidity about S0% at 
220 K) would appear as a feature about 10°K colder than the surrounding scene. 

4.3.2. Observed effects 

UW has been monitoring VAS SST along the East Coast during 1982. They report a few degree 
depression in the VAS SST at buoys located east and west of the Florida peninsula during June 1982, 
and smaller depressions are observed during mid-summer which are consistent with AVHRR SST 
depressions. 

Because VAS observations are oriented towards the United States, they rarely extend below 
30°N. By chance, some 3-channel low-latitude scans are available from the VAS Demonstration at 
GSFC, taken at 1530, 1600 and 1630 GMT on 6 and 8 April 1982, which is a few days after the 
main eruption of El Chichon. 

a) On both days, an unusual low cloud was observed in the 1 1 and 1 2 pm VAS radiances. 

The cloud was about 10°K colder than the surrounding clear areas, and it obscured the 
quadrant northwest of El Chichon. Fig. 4.8 shows this cloud’s effect on the VAS SST 
calculated over the entire scene from the split window radiances observed at 1 600 GMT 
on 6 April 1982. El Chichon’s location is marked by a “+”. The “SST” is 5°C less than 
the surroun >lngs near Mexico and over the Gulf of Mexico. There is no clear sign of an 
effect on VAS ‘'SST” further west of El Chichon, where stratospheric particulates might 
be expected. 

b) Fig. 4.9 shows the corresponding 6.7 pm image. There is an unusual cool (7°K lower) 
feature in the warm cloud-free streak over Mexico, to the west of El Chichon. This molec- 
ular water vapor feature is more compact than the contemporaneous stratospheric features 
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Figure 4.8. The sea surface temperature measured by VAS from GOES-5 at 
1600 GMT on 6 April 1982. The SST grey scale is at the top. 
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observed in the SBUV/TOMS images (Section 4.4). One millimeter of precipitable water 
over an area of 2S0 by 400 km would be 10 s tons of HjO. The feature is net identifiable 
in the 6.7 pm image on 8 Ap^"!. 

Without a more continuous series o r VAS images, one cannot be certain that either of the VAS- 
observed features was actually ejected ’ey El Chichon. 

4.3.3. Corrections 

Beca ise the VAS instrument is normally dedicated to studies over the United States, where 
the effects of El Chichon are comp^rajle to VAS’ radiometric errors (±1°K brightness temperature), 
there is no need to correct VA C 'ings for El Chichdn’s ejecta. 

4.3.4. Further research 

Future volcanic erruptions in the Western Hemisphere should be monitored more deliberately 
with short (10 minute) interruptions of the operational V1SSR schedule to acquire 12-channel VAS 
soundir js over a 10° latitude band every hour. A volcano watch of this sort would be basically the 
same as hurricane and severe storm watches already permitted by NOAA. During such a volcano 
watch, the VAS Demonstration facility rt GSFC should collect VAS radians and provide immedi- 
ate data processing for scientific investigators. 


D. Clusters 
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4.4. TOMS AND SBUV TOTAL OZONE RETRIEVALS AND VERTICAL PROFILES 

4.4.1. Total ozone determination 

4. 4.1.1. Theoretical effects 
a. Gaseous absorption 

The TOMS and SBUV instruments measure total ozone by observing the radiance of the back- 
scattered atmospheric sunlight in the Hartley absorption bands of ozone. In the unperturbed 
atmosphere other known species which could interfere with total ozone retrievals have abundances 
<1(T 4 of the ozone abundance and do not cause retrieval errors Volcanic clouds contain two 
sulfur compounds, SO2 and CS2 , which can directly interfere with ozone if the quantity is great enough. 
The spectra of SO2 and CS2 from Wu and Judge ( 1 98 1 ) are shown in comparison with ozone in Fig. 4-1 0. 
Three of the wavelengths used by SBUV and TOMS for total ozone measurement are indicated by the 
arrows. The retrievals use the wavelength pairs 3 1 2.5/33 1 .2 nm and 3 1 7.5/339.8 nm, reierred to as the A- 
and B-pairs, respectively. It is clear that either constituent will produce erroneously high ozone 
amounts if not corrected. The absorption coefficients (base e) for O3, S0 2 and CS 2 at the total 
07^>ne wavelengths are listed in Table 4.4 together with the ratios of S0 2 and CS 2 to ozone absorp- 
tion coefficients. The differences in the absorption coefficients for the A- and B-pairs appear on the 
last two rows. Although S0 2 and CS 2 are expected to have significant temperature dependencies, 
measurerr nts at low temperatures corresponding to El Chichon cloud exist only for 
(Wine et al., 1981). The effects of S0 2 are seen to be strong at the shortest two wavelengths which 
are the primary ones for the ozone measurement. One unit of S0 2 is equivalent to about 2.5 units 
of ozone. A unit of carbon disulphide is roughly equivalent to that of ozone. At 331.2 and 339.8 
nm S0 2 absorption is small compared with ozone absorption while CS 2 is considerably larger. 

The true S0 2 cross sections at El Chichon cloud altitude of 27 km will be smaller than the above 
values due to the temperature effects, but the amount is at present unknown. 
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Figure 4.10. Cross section of ozone, sulfur dioxide and carbor. disulphide in the 
spectral region used by SBUV and TOMS for total ozone measurements. 
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Table 4.4 

Absorption Coefficients (Base e) at TOMS and SBUV Wavelengths 
in Units of cm -1 of Gas at NTP 


Wavelength 

(nm) 

a O3 
(229 K) 

0 SO2 
(294 K) 

aCS 2 
(254 K) 

a SC>2 
a0 3 

a CS2 

a O3 

312.5 

1.777 

4.09 

1.602 

2.30 

0.902 

317.5 

0.9264 

2.36 

1.274 

2.55 

1.375 

331.2 

0.1676 

0.042 

0.632 

0.251 

3.771 

339.8 

0.0405 

0.019 

0.080 

0.469 

1.975 

360.0 

0 

0 

0 

- 

— 

380.0 

0 

0 

0 

- 

- 

A-pair 

1.609 

4.048 

0.970 

2.52 

0.603 

B-pair 

0.8859 

2.341 

1.194 

2.64 

1.348 


b. Aerosol scattering 

The present retrieval model assumes a multiple scattering Rayleigh atmosphere with a 
Lambertian reflecting lower boundary. The peak contribution to backscattered radiances is found 
in the lower troposphere at tropical latitudes. A dense aerosol cloud in *he stratosphere can scatter 
photons to space without having passed through the o^one below the cloud and thus produce an 
underestimate of the total ozone. However, u the particles also absorb, this effect is counteracted. 
Dave (1978) has studied these effects for typical stratospheric and tropospheric aerosol distribu- 
tions. He finds very slight ozone underestimations for pon-absorbing aerosols but significant over- 
estimations for absorbing aerosols. Since El Chichon cloud is higher and denser than Dave’s 
assumed aerosol distributions, it will be necessary to carry out radiative transfer calculations to 
determine the particular effects. 

4.4.1. 2. Observed effects 

El Chichon cloud produced a strong perturbation in T OMS radiar.ces at the two shortest 
wavelengths (Krueger. 1983). This is shown in Fig. 4 11, which portrays the logarithm cf me ratio 
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Figure 4. 1 1 . Apparent total ozone amounts in a trace across El Chichon cloud (upper). 
Equivalent optical depth of atmosphere in same trace across cloud at five TOMS wave- 
lengths (lower). 
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of radiance (I) to solar irradiance (F) at five of the TOMS wavelengths for a trace across the early 
cloud. The spectral signature is consistent with the presence of S0 2 and inconsistent with CS 2 , 
as can be seen by inspection of Table 4.4. The longer wavelengths show only slight perturbations 
due to aerosols. Thus by far the largest effect is a severe overestimation of total ozone by more 
than 200 m atm cm due to the S0 2 absorption (see Fig. 2.6 in section 2). As the cloud was sheared 
by the winds and the S0 2 was chemically lost, the ozone retrieval error diminished. 

A second effect was that the difference between A- and B-pair ozone retrievals was large over 
the cloud. If the S0 2 absorption coefficients in Table 4.4 were correct and the dust was non- 
absorbing, the difference should be small. An explanation will require low temperature S0 2 cross 
section data and radiative transfer calculations. 

A third effect, found after about 2 months was scan angle dependence of the TOMS ozone 
retrievals with lower values obtained at nadir than at large scan angles. This effect is not understood 
but probably is due to aerosol effects. 

4.4. 1.3. Correction of data 

The total ozone retrieval errors due to the presence of S0 2 diminish rapidly as the volcanic 
cloud is dispersed by winds. An upper limit to the total ozone error is obtained by assuming no 
chemical conversion, while the cloud is mixed uniformly within a latitude zone. For example, if 
El Chichon injection of S0 2 is mixed uniformly within the latitude zone 0°-30°N, the average bias 
in the estimate of ozone is 2.3 m atm cm. Given a determination of the rate of loss of SO->, it 
should be possible to cotrect ihe zonal mean total ozone with adequate precision for most purposes. 

4.4.1 4. Further research 
a. Analysis of existing data 

The properties of the cloud will evolve with time. Effort is needeJ to determine the best cloud 
signature and to use this to locate both the tropospheric and stratospheric clouds as a function of 

4-34 


-) 



time. This will permit a determination of the rate of loss of S0 2 which is required for any correc- 
tion scheme. Aerosol properties may also be inferred by comparison with surrounding areas during 
the time when a discrete cloud can be identified. 

b. Comparison with other satellite data 

The cloud boundaries differ between those derived from TOMS, AVHRR, and VAS observa- 
tions. This may be due to wavelength effects but also could result from altitude differences across 
the cloud. This needs to be explored as a possible way of discriminating cloud altitude. In addition 
the SBUV continuous scan data need to be examined for similar reasons as noted in the ozone 
profile section and also to determine whether constituents other than S0 2 were present. 

c. C oss section data 

Sulfur dioxide cross sections at stratospheric temperature and pressure conditions need to be 
men ured. 

d. Radiative transfer calculations 

The conditions of an absorbing-scattering aerosol -loud near the ozone maximum have not 
been simulated. Off-nadir radiative transfer calculations have not been made for TOMS under high 
aerosol conditions. These calculations are essential for understanding El Chichon effects. 

A. J. Krueger 


4.4.2. SBUV ozone profiling 

4.4.2. 1. Theoretically expected effects 

The effect is very sensitive to the altitude of the volcanic cloud - the higher the cloud the 
larger the effect. Aerosols that rise above the ozone density peak (~ 25 km in tropics) produce a 
far greater effect than those that lie below the peak. Therefore, the effect of El Chichon on SBUV 
profiling should greatly exceed that of Mt. St. Helens, which injected aerosols in the lower 
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stratosphere only. If El Chichon cloud was located at 27 km, the four shortest SBUV bands, 
located below 290 nm, will not be affected. The effect at 292 nm is strongly dependent on the 
aerosol altitude profif and the distribution of ozone in or above the aerosol layer — the higher the 
aerosol layer in relation to the ozone layer the larger the effect. If we assume that little aerosol 
was injected above 30 km, the effect on 292 nm radiances should be less than 5% (assuming trop- 
ical ozone profiles). For longer SBUV profiling bands located at 298, 302 and 306 nm the predicted 
effect of the aerosols is extremely large ~ 30-50%. If the SBUV profiling algorithm ignores the 
presence of such a large perturbation, the ozone profiles derived in 1-5 mb region would be in error 
by 5-10%; below 5 mb the error would increase rapidly, perhaps becoming as large as 50% in the 
10-20 mb region of the atmosphere. These results are, however, preliminary and more 
work needs to be done to understand the effects. 

Dave (1978) made an extensive modeling study of the effects of stratospheric aerosols on the 
estimates of total amount of ozone as measured by SBUV. His aerosol model, however, was not 
the same as El Chichon aerosol. Dave’s aerosol profile had a maximum at 20 km, whereas 
El Chichon ’s maximum concentration was at 26 km. Also, Dave’s maximum aerosol optical thick- 
ness was 0.06. as compared with an El Chichon value of 0.2 (Fig. 3.2). This study indicated that 
for non-absorbing aerosols the effect on estimates of total amount of ozone is small but can become 
quite significant for absorbing aerosols. 

The s .ape of the scatterir i <ase function strongly affects the amount of scattered light, 
but it is fairly well understood. Since the SBUV is a nadir viewing instrument, the scattering angle 
is the same as 180° minus the solar zenith angle. Based on calculations done by M. King and 
John DeLuisi (private communication), the largest effect would be for overhead sun and the mini- 
mum e f fect should occur at a solar zenith angle of ~50°, with variations of the order of a factor of 
five. Some interesting anomalies have been predicted for the aerosol phase function - showing 
oscillations at scattering angles of ~ 160° - that should be observable in SBUV radiances. 
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The single scattering albedo of the volcanic aerosols for the ultraviolet spectrum is, perhaps, the 
greatest unknown in assessing the impact of the volcano on SBUV and TOMS. Even a weak absorption, 
if present, would significantly reduce the effect of aerosol on SBUV radiances for the shorter wavelengths. 
Indeed, for an aerosol with imaginary part of index of refraction of ~0.05, the aerosol might appear 
neutral at longer wavelengths-producing no observable signatures in the back-scattered radiances. 

The presence of SO-, will also cancel out the effect of aerosols and could, indeed, reverse the 
effect - producing a decrease of radiances. However, preliminary results indicate that the strato- 
spheric SO-, is not large enough to cancel out the scattering effects. Even if it did at some wave- 
length, the effect would vary dramatically with wavelength, following the unique SO-, absorption 
spectrum. 

4.4.2. 2. Further research 

No study has so far been done that predicts the effect of an aerosol layer located above the 
ozone density maximum on the SBUV radiances - particularly those that are strongly multiply scat- 
tered. Existing studies were either done for a very thin aerosol layer at ~20 km (Dave) or for clouds 
located in the troposphere (Dave. Fraser). These studies cannot be used to predict the effect when 
the aerosol layer is above the ozone layer. 

The effect of a large perturbation of the longer wavelength radiances on the upper level profile 
retrieval is not known and should be studied. 

It is extremely impoitant for the analysis of SBUV/TOMS data that we know the single scat- 
tering albedo of El Chichon aerosols on ultraviolet wavelengths. Some measurements already exist 
but are not sufficient (Patterson. 1981). 

Optical properties of the volcanic aerosol can be obtained by analysing the SBUV/TOMS data. 
Such a study could reveal the phase function of the aerosols, its optical depth, altitude distribution. 
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and amount of SC >2 present. More importantly, if one can find a way to track the cloud over the 
globe and the way is disperses, it may be possible to correct the SBUV and UMKEHR ozone data 
using those results. 

P. K. Bhartia 
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4.5. ERB ON NIMBUS-7 

The Earth Radiation Budget (ERB) experiment on Nimbus-7 measures upward flux in two 
broad wavelength bands: 

1) 0.2 - 50 pm 

2) 0.2 - 3.8 pm 

The instantaneous fie'd-of-view extends beyond the earth’s limb. Thp difference between total 
and solar upward fluxes i- the long wave upward flux. The major influence of stratospheric 
aerosols is on the solar reflectivity, which is anticipated to increase slightly for regions covered by 
El Chichon aerosol cloud. The total change in reflected energy is expected to be 5-10 W/m 2 for 
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El Chichon compared with 200-250 W'm 2 mean solar reflection, and this enhancement could pos- 
sibly be detected if the time series is analyzed carefully. Unfortunately, ERB data processing is 
currently at 1980, and data pertinent to El Chichon will not be available for at least another year. 
There is also no gain to be expected from analyzing the data out of sequence, since comparison with 
previous data is the only way to detect the order of the change anticipated. 

The effect of the cloud in the long wave is more difficult to assess. However, since the optical 
depth of H 2 S0 4 in the long wave is much less than in the solar and the high aerosol cloud is at a 
temperature close to the middle troposphere, no significant effect is anticipated. 

Figure 4. 1 2 shows the additional reflection of solar energy to bt expected from a layer of 157c 
H 2 S0 4 aerosols of visible (0.55 pm) optical depth of 0. 1 . For small optical depths, the response 
is almost linear and the figure can be used to scale the effect of El Chichon cloud, which had an 
optical depth of 0.2-0.25 in a band extending from 0°N-30°N in the summer of 1982. 

The cloud is expected to diffuse over time and the radiative perturbation is expected to 
decrease. It must be stressed that the Fig. 4. 1 2 only supplies the estimated initial radiative pertur- 
bation. If through some physical response, there is a significant change in tropospheric cloudiness 
as a result of this radiative j. rturbation, El Chichon signal will probably not be detected in the 
zonal mean reflectivity. 

Although the radiative perturbation is only marginally detectable, the aerosol cloud could 
significantly influence the mean surface temperature. Application of the radiative model to the 
GLAS Multi-Layer Energy Balance Model with an assumed two-year residence time for the 
aerosol, indicates that the peak response of climate should occur approximately 2.5 years 
after the aerosol layer has reached its production peak. The peak effect on the hemisr ric 
mean surface temperature will be between 0.3 and 0.5°C depending on the assume parameters 
and climate feedback processes. The low latitudes will sho'v the earliest response if the aerosol 
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layer is largely confined to those latitudes, as now seems to be the case. The.e will then be a 
cascading effect on higher latitudes, with the peak response appropriately delayed. 

Hitshvardhan 
M. D. Chou 
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4.6. CZCS PRODUCTS 

The data processing for the Coastal Zone Color Scanner (CZCS), as part of the «oi.eme to 
eliminate atmospheric path radiance consisting of both Rayleigh and aerosol contributions, begins 
by calculating the Rayleigh radiance for the bands at 443, 520, 550 and 670 nm, and subtracting 
that radiance from the aperture radiance measure-’, by the sensor. The algorithm for the atmospheric 
conection then utilizes an assumed aerosol scattering efficiency of 1 for the four previously men- 
tioned wavelengths, and the assumption tliat the reflected radiance from the surface at 670 na. 
meters is negligible in the case of clear water. The algorithm searches for the clearest or bluest water 
and then calculates a relative scattering efficiency of the aerosol at ^43, 520 and 550 based on a 
relative scattering efficiency of 1 at 670 nanometers. 

Although aerosol path radiants arc calculated for all four spectral bands in the removal of the 
atmospheric contributions of CZCS data, only the aerosol path radiance at 670 nanometers is 
retained and archived as part of the Level II r>T oduct. The aerosol path radiance for ti.e channel r ' 
670 nanometers is retained (as part of the pictorial and the magnetic tape product) in a range of 
0 to 1.5 milliwatts per square centimeter micrometer steradian. This path radiance has shown that 
aerosols are quite spatially variable and that extrapolation from one area of the picture to another, 
such a- i, - i water to adjacent land, would probably not be valid 
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At the southern end of Baja California in the latitude range most affected by the El Chichon 
dust, an aiea has been selected and a number of scenes have been requested for priority processing 
through to Level II, so that the aerosol path radiance in this area might be compared wii from 
previous years in the same area. Since tne CZCS on Nimbus-7 w:.s launched in 1978, we have a 
data set from the summer of 1979, 80 and 81, prior to the El Chichon eruption, that can provide 
a comparison of the aerosol path radiance for three years running to that aPer the El Chichon 
eruption. The data from the orbits in question is being processed now as rapid, , „ possible by the 
Information Processing Division at NaSA/GSFC, and the comparison will be made as scon as the 
Level II product is available. 


W. A. Hovis 


4." SMMR PRODUCTS 

The Scanning Multichannel Microwave Radiometer (SMMR) on Nimbus-7 is a ten-channel 
microwave radiometer that measures brightness temperatures at 6.6. 10.7, 18, 21, and 37 GH7 in 
dual polarization. The conical scanning geometry of SMMR is such that the i strument views the 
earth’s surface at a constant incidence angle of about 50°. The spatial resolution of the instrument 
is proportional to wavelength and varies from about 1^5 .:m a, 6.6 GHz to about 30 km at 37 GHz 

Except possibly for a eriod within 1 7 hours of a major El Chichon eruption, ejecta will not be 
discernible by satellite microwave observations. To explain the negligible efl’Ct of the El Chichon 
aerosols, a comparison can be made with the empirical insults of SMMR observations of liquid 
waier clouds. The microv ave index of retraction of water is large, and one anticipates the index 
o f refraction of El Chichon aerosols that reside in the stratosphere a few months after the last 
eruption is lar.'e, also. The cloud drops are several orders of magnitude smaller than the -microv ave 
wavelength (~ l cm), and El Chichon particles are even an or^e: of magnitude smaller than the cloud 
drops (section 3.2). Attenuation of microwave radiation by suen clouds is dominated by absorption. 
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rather than by scattering (van de Hulst, 1957). The total absorption is proportional to the cloud 
mass in a vertical column. The 37 GHz, horizontal polarization channel of the SMMR is the most 
sensitive to clouds over water. Since the noise-equivalent-temperature of the channel is about 1 K, the 
minimum detectable cloud mass is about 1 mg cm -2 (Chang and Wilheit, 1979). The mass of El 
Chichon cln»» ’ a few months after the last eruption was 0.006 mg cm -2 (section 3.2), which is far 
too small to oe detected by SMMR. 

The sea surface temperature derived from the SMMR observations can possibly be used to 
verify the NOAA operational SST derived from infrared radiances. The SMMR SST can be obtained 
globally for monthly periods. Then, the deviation between the microwave and infrared SST*s can 
be used to outline the geographical extent of the El Chichon stratospheric cloud, since the infrared 
radiances are affected by the cloud, whereas the microwave radiances are not. 

P. Gloersen 
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4.8. MSS AND TM ON LANDSAT 

The Multispectral Scanner (MSS) and the Thematic Mapper (TM) measure the radiance of 
sunlight scattered from the earth and its atmosphere in many spectral bands between 450 to 2300 
nm. Characteristics of surface features on the earth are derived from these observations. The major 
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cause for atmospheric error in classifying surface features is the variable optical characteristics of 
atmospheric aerosols. An increase in amount of aerosols reduces contrast, reduces spatial resolution, 
and alters spectral signatures. 

The effect of El Chi chon aerosol on MSS and TM observations is expected to be less than 
effects due to the atmosi>here that are routinely encountered. The aerosol effect can be related to 
its optical thickness. A few months after the last major eruption. El Chichon’s stratospheric aerosol 
had stabilized with an optical dcp’h of 0.2 for the spectrum 400-1000 nm in the latitude zone cf 
10° -30°N. If the aerosol diffuses over the entire northern hemisphere, the optical depth will be 
0. 1 , or less. Its gradients would be weak with an expected uniform spatial distribution. 

The dominant source of atmospheric error in characterizing surface features, however, would 
still be the tropospheric aerosol, concentrated within a few kilometers of the earth’s surface. The 
mean value of that aerosol optical thickness is 0.5 and its standard deviation is about one-half of 
that over the eastern half of the United States during the summer. The stratospheric cloud would 
add an additional thin, uniform veil between Landsat and the ground. This veil is not expected to 
cause additional effects that are important in comparison with tropospheric effects that occur 
routinely. 



5. CONCLUSIONS 

The last eruption of the El Chichon volcano occured on 4 April 1982 (day 94). The AVHRR 
and VISSR delineated the diffusion of the aerosols resulting from the volcano as they circled the 
globe during the next 3 weeks. The AVHRR sea surface temperature bias indicated that the strato- 
spheric aerosols were mostly confined to the latitude band 10N-30N before diffusing north and 
south in July and later. SME and lidar observations placed the aerosol layer between 20 and 35 km 
altitude, with a weak layer between 1 8 and 20 km. 

The effect of the stratospheric cloud on measurements of sea surface temperature depends 
on the sensor and its algorithm. Temperatures derived from the AVHRR observations were as much 
as 2.5°C lower than ship and buoy temperatures. The sea surface temperatures have not yet been 
derived for other infrared measurements (H1RS and VAS). Because the stratospheric cloud is es- 
sentially transparent to microwave radiation, SMMR and MSU measurements will not be affected 
by the cloud. 

TOMS observations of the total amount of ozone in the region of the volcanic ejecta were 
seriously over estimated because of the absorption by SO-,. Krueger, however, has utilized the 
TOMS observations to estimate the total amount of SO-, ejected by El Chichon. Both the strato- 
spheric aerosols and S0 2 are expected to have significant effects on the vertical profiles of the 
ozone concentration, derived from SBUV observations. 

The temperature and humidity profiles measured by the infrared sensors will be affected by 
the volcanic ejecta, at least in the more concentrated regions. VAS profiles have been restricted to 
the United States. The profiles already analyzed were taken before the cloud reached the U.S. 

When the cioud reaches the U.S.. it would be less concentrated. Then the temperature perturbations 
are expected to be less than 1 K, which is within the noise level. Chesters feels that the VAS tem- 
perature and humidity retrievals over the United States will not be affected by the stratospheric 
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cloud. Susskind, on the other hand, has seen large positive and negative temperature anomalies in 
the H1RS observations within a day of an El Chichon eruption and near to it. Small biases in the 
temperatures that persist after the stratospheric cloud becomes more tenuous may still have a 
significant effect on the retrievals. 

The stratospheric cloud has a negligible effect on microwave observations and is expected to 
oe small or routinely correctable for ERB, CZCS, MSS, and TM satellite observations. 

In order to modify the satellite algorithms to account for th v volcanic effects, additional work 
is required. The anomalies can be used, however, to measure the characteristics of the volcanic 
aerosols and dispersion. 
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6. RECOMMENDATIONS 

Algorithms for some of the satellite sensois need to be modified in order to account for the 
effects of volcanic ejecta on radiant energy. The SBUV/TOMS measurements in the ultraviolet 
spectrum of the total ozone and its vertical profile are strongly affected by the residual stratospheric 
aerosols and trace gases. The radiation characteristics of the aerosols in ultraviolet light are 
poorly known, and the existing SO2 absorption spectra are not accurate. Experimental data on both 
scatterii g and absorption characteristics are needed. With such data the radiative effects of strato- 
spheric clouds need to be modeled. Then the effects of the aerosols can be accounted for in not 
only the ozone algorithms, but also for the infrared algorithms. It is recommended that the above 
studies and experiments be carried out. 

The residual cloud remaining in the stratosphere causes an apparent decrease in the sea surface 
temperature measured in the infrared spectrum. Stowe is developing an algorithm to remove such 
a bias from the AVHRR data. Removal of this bias from the H1RS data presents a different problem 
and needs to be analyzed. SMMR measurements of SST are not affected by the volcanic cloud and 
can be used as a set of reference data. It is recommended that studies be made to account for 
El Chichon cloud on SST. 

The temperature and humidity sounding algorithms for HIRS measurements require modifi- 
cation to account for the volcanic ejecta. Susskind intends to analyze the HIRS data carefully. 

His studies will indicate the requirements to derive accurate temperature and humidity profiles. 

It is recommended that such studies be made, and Susskind has indicated his intent to do so. 

The satellite data already available can be used to provide more information about the effects 
of El Chichon eruption on the atmosphere. The amounts of trace gases and their vertical profiles 
may be derived from the SBUV measurements of the radiance with high spectral resolution in the 
spectrum X 1 600-4000A. Krueger has already made preliminary estimates of the amount of SO-, 
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injected into the atmosphere from the TOMS data. Both SBUV and TOMS may possibly provide 
aerosol optical thickness data about aerosol particles smaller than a few tenths of a micrometer in 
diameter. Since volcanic aerosol effects have not been apparent in the TOMS data, all radiation 
parameters will have to be carefully accounted for in order to measure the aerosol effects. It is 
recommended that research be performed on the information that can be derived from satellite ultra- 
violet observations concerning volcanic ejecta, its composition, spatial extent, and decay characteristics. 

The visible and near infrared radiances measured by AVHRR and VAS over water can be used 
to measure the geographical dispersion of aerosols. Stowe is investigating the possibility of doing 
this. Also, the geographical location of the stratospheric aerosol layer closely coincided with the 
region of sea surface temperature anomalies measured by AVHRR. The length of time such anom- 
alies exist and their relation to the stratospheric layer are under investigation by NOAA scientists. 

It is recommended that such studies be continued concerning use of visible and near infrared 
observations to characterize the volcanic aerosols. 

The satellite data, together with additional data from surface, aircraft, and balloon observations, 
provide a rich source to utilize in investigation of climatic effects, atmospheric diffusion, and the 
complex chemistry of the stratosphere. 


6-2 



